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We present XMM-Newton Reflection Grating Spectrometer observations of X-ray 
clusters and groups of galaxies. We demonstrate the failure of the standard cooling- 
flow model to describe the soft X-ray spectrum of clusters of galaxies. We also em- 
phasize several new developments in the study of the soft X-ray spectrum of cooling 
flows. Although there is some uncertainty in the expected mass deposition rate for 
any individual cluster, we show that high resolution RGS spectra robustly demon- 
strate that the expected line emission from the isobaric cooling-flow model is absent 
below 1/3 of the background temperature rather than below a fixed temperature in 
all clusters. Furthermore, we demonstrate that the best-resolved cluster spectra are 
inconsistent with the predicted shape of the differential luminosity distribution and 
the measured distribution is tilted to higher temperatures. These observations create 
several fine-tuning challenges for current theoretical explanations for the soft X-ray 
cooling-flow problem. Comparisons between these observations and other X-ray mea- 
surements are discussed. 



1. The Cooling-Flow Model 

A long-standing theoretical prediction is that the 
intracluster medium in the cores of galaxy clusters 
should c ool by emitting X-ray s in less than a Hub- 
ble time (iFabian fe Nulsen l ll977t ICowie fe Binnevlll977l 
iMathews fc Breg manlll9 78l). If the cooling proceeds in- 
homogenously ijNulserJ |l986t iJohnstone et al.lll992|) . a 
range of temperatures is likely to exist in the cool- 
ing plasma. Straightforward thermodynamic arguments 
show that at constant pressure the particular distribution 
of plasma temperatures is described by the differential 
luminosity distribution, 



dL _ 5 M 
dT 2 [im v 



(1) 



where k is Boltzmann's constant, [im v is the mean mass 
per particle, and M is the mass deposition rate. This 
yields a unique X-ray spectrum for an assumed set of el- 
emental abundances and a given maximum temperature 
where the cooling is assumed to begin. Inhomogenous 
cooling of this sort is broadly called a cooling flow. 

We can apply Equation 1 to the X-ray spectrum 
by using standard plasma codes to predict the energy- 
dependent luminosity at a given temperature. There are 
three major assumptions: First, it is normally assumed 
that the total luminosity emitted from the high temper- 
ature plasma is not substantially different than that ob- 
served in X-rays. If this were not true or if the plasma 
was indirectly transfering its thermal energy to another 
source, then this would invalidate the model. Second, we 
assume that there is no substantial heating that would 



modify the thermodynamic history of the plasma. Third, 
we assume that the bulk of the cooling volume is ap- 
proximately in a steady-state, so that the cooling flow is 
not dramatically changing during a characteristic cool- 
ing time. Although many modifications and theoretical 
challenges to this simple model have been discussed for 
more than 25 years, it was not known until recently if 
the X-ray spectrum deviated significantly from this sim- 
ple cooling-flow model. 

In order to measure a differential luminosity distribu- 
tion at keV temperatures like that in Equation 1, emis- 
sion lines from individual ions need to be measured. In 
particular, it is important to measure emission lines from 
the Fe L series, which are ions that happen to have their 
ionization potentials in the temperature range near and 
below cluster virial temperatures. This is shown in Fig- 
ure 1 where the relative charge state abundance (ele- 
mental abundance times fractional ionic abundance) of 
various ions is plotted. However, measuring the relative 
strength of emission from individual Fe L ions is only pos- 
sible to do in detail with resolving powers (^j) above 
100. These resolving powers are not achievable with 
non-dispersive CCD experiments like ASCA/SIS, XMM- 
Newton/EPIC, and Chandra/ACIS, although some in- 
formation can be obtained from the shape of the unre- 
solved Fe L emission line complex. The Reflection Grat- 
ing Spectrometer (RGS) on XMM-Newton, however, is 
capable of high-resolution spectroscopy for sources with 
spatial extent similar to many cooling flows. 

The RGS is a dispersive spectrometer that contains 
a grating array that is placed behind the XMM mirrors. 
The grating array deflects photons to a CCD bench. The 



2 



Peterson et al. 



ur' J 










io- 4 








S XVI & XV : 

Si xra & xiv : 

Mg XII & XI 

Ne X & IX 

ovm&vn 










N VII & VI : 
C VI&V 1 




: /7 \ \ \ 








in*" 8 


/ \ A 






1 


to- 8 


O.l 


1.0 


FeXVU 
FeXVI 
FeXV 


FeXXVI - 

Fe XXV = 

FeXXTV 

Fe XXIII '. 


IO"* 






FeXIV 
Fe Xni 
FeXn 
FeXI 
- FTX 


FeXX ■ 
"TeXCt — ■ 

Fe XVin 


10" S 






Fe IX 





Temperature (keV) 

Fig. 1. — Fractional abundance of a given ion plotted 
again st temperature in keV (taken from lAr naud & R aymond! 
1992). The fractional abundance is multiplied by the abundance 
of that element relative t o hydrogen in the solar neighborhood 
lAnders fe G revessc 1989). The top plot shows the helium and 
hydrogen-like charge states for several elements. The bottom plot 
shows a number of charge states of iron. 



RGS was designed to have a relatively high dispersion 
of 3 degrees for the soft X-ray band in order to com- 
pensate for the fact that the X MM mirrors would blu r 
sources by about 10 arcseconds l|Kahn fe Hettrickl ll985'). 
Cooling flows happen to have a similar size to the XMM 
blurring, so that observations of cooling flows also ben- 
efit from the large dispersion angles. The dispersion of 
an incident X-ray is determined by both the wavelength 
of the photon and the angle it hits the grating array. 
Therefore, for a source with significant spatial extent 
like a cluster, the coupling of the sky position with the 
wavelength has to be modeled. A Monte Carlo approach 
has therefore been used to analyze extended source RGS 
observations. Some discuss ion of this can be found in 
iPeterson. .Ternigan fc Kahnl I 



2. Spectroscopic Constraints on Cooling Flows 

The initial application of the standard cooling-flow 
model to a high resolution spectrum of a cluster of 
galaxies is shown in Figure 2. The top panel shows 
the spectral prediction for a model based on Equation 
1. The second panel shows the comparison between the 
model (blue) and the spectrum of Abell 1835 (red). No- 
tably absent in the data are Fe XVII lines. Finally, 
the bottom panel shows a model where the emission 
below 3 keV is suppressed. This fairly simple model 
of the X-ray emission seems to describe the data quite 
well. The plasma appears to match the cooling flow 
model betwen 3 keV and the maximum cluster tem- 
perature of 8 keV but not below 3 keV. The obvious 
theoretical question is why it appears close to match- 
ing the cooling flow model at high temperatures, but 
not at low temperatures. It is clearly not an isother- 
mal spectrum, but it is also inconsistent with a complete 
cooling flow model. Full de t ails of this comparison can 
be found in IPeterson et al.l l)20f)l[). Other grating ob- 
servations contain similar results llTamura et al.l l2001e 
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Fig. 2. — Comparison of the cooling flow model to Abell 1835. 
The top panel shows the cooling flow model based on Equation 
1 in the text. Note the numerous Fe L lines predicted in the 10 
to 17 A region. The middle panel shows the comparison of that 
model with the data from the galaxy cluster Abell 1835 which was 
thought to host a massive cooling flow. The bottom panel shows 
the cooling flow model with the emission suppressed below 3 keV. 
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Fig. 3. — The emission spectrum of the standard cooling flow 
model broken up into several temperature bins. The red line shows 
the 3 to 6 keV spectrum, the yellow shows the 1.5 to 3 keV spec- 
trum, the green shows the 0.75 to 1.5 keV spectrum and the blue 
shows the 0.375 to 0.75 keV spectrum. The normalization of each 
of the components is adjusted to see where the cooling flow model 
fails. 



If we decompose the spectrum in the top panel in Fig- 
ure 2 into different temperature ranges, we arrive at the 
spectra in Figure 3. Here we plot the 6 to 3 keV (red), 3 
to 1.5 keV (yellow), 1.5 to 0.75 keV (green), and 0.75 to 
0.375 keV spectrum (blue). Then if we fit the observed 
X-ray spectra while allowing the relative normalization 
of each of these four spectra to vary, we can fit the X-ray 
spectra of 14 cooling-flow clusters of galaxies quite well. 
This is shown in the spectral fits in Figures 4-6. These 
fits also use a spatial model for the cluster emission folded 
through the Monte Carlo for the RGS. The red line shows 
the fit, the blue histogram is the data, and the green line 
is the standard cooling flow model with no adjustments. 
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Fig. 4. — Fluxed spectra of several hot (about 5 keV) clusters of 
galaxies. The red line is an empirical fit allow the normalization 
of the cooling flow model for several temperature ranges to be 
adjusted. The green line is the standard cooling flow model with 
no adjustments. The model fails most strongly from Fe XVII lines, 
which implies that the lowest temperature emission is absent. 



The cooling flow model deviates significantly from the 
data in a similar manner to the Abell 1835 comparison 
in Figure 2. The low temperature emission is not present 
or significantly reduced in the cluster spectra. Full de- 
tails of the model and analysis method can be found in 
iPeterson etaD l|2003bft . 

Figure 7 shows the differential luminosity of each of 
the four spectra (proportional to the total normalization 
given by the fit) against the temperature. Each cluster 
is shown in a different color and the actual detections 
are connected by a line. If the cooling flow model were 
correct, then the points should line close to the line y — 1. 
Figure 8 shows the same plot but it is plotted against the 
temperature relative to the maximum temperature. The 
points appear roughly consistent with a power law in 
fractional temperature where the index is near 1 to 2. 

Figure 7 and 8 demonstrate several aspects of the fail- 
ure of the cooling-flow model. First, it is clear that sig- 
nificant quantities of plasma below the maximum tem- 
perature always exists. Second, the model appears to fail 
at a fraction of the maximum temperature rather than a 
fixed value as shown by the luminosity distribution cross- 
ing the line y — 1 near T = | in Fi gure 8. Furthermore, 
the model fails in its shape (it is not horizontal) since it 
is tilted towards higher temperatures and the normaliza- 
tion at low temperatures is not the only discrepancy. The 
overall normalization, however, is somewhat difficult to 
interpret in the absence of a new theoretical model, since 
a new model would probably add or subtract heat to the 
system which affects the normalization. There is also 
significant scatter in both plots. It is not clear whether 
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FlG. 5. — The same plot as Figure 4 but with medium tempera- 
ture (3 to 5 keV) clusters. 



l|i , Abell 262 


!; ib is 


20 25 30 35 — 

Abell 1837 


io i l 




20 25 30 35 

M87 


I 10 1 




20 25 30 35 

NGC 533 _ 


10 . 15 


20 25 30 35 

Composite Spectrin] 



10 15 20 25 30 35 

Rest Wavelength (A) 



FlG. 6. — The same plot as Figure 4 and 5 but with cool (less 
than 3 keV) clusters and groups. 



this is due to subtle fitting degeneracies or is a real dif- 
ference between cooling flows. We also cannot be certain 
if the empirical model, the power law in fractional tem- 
perature, continues to arbitrarily low temperatures. It is 
clear, however, that the model is off by as much as an or- 
der of magnitude at the lowest temperatures. Finally, it 
appears that the simple temperature cut-off model used 
in Figure 2 might be oversimplified. The best-resolved 
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Fig. 7. — A comparison of the luminosity in specific temperature 
ranges to the cooling flow model. See the text for details. 
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Fig. 8. — Same as Figure 7 but plotted as a function of the 
fraction of the maximum temperature. 



spectra like that for 2A0335+096 has emission lines from 
Fe XXIV through XVIII. This has forced us to use a more 
complex model than the cut-off model to match the low- 
est temperature emission. Similarly, just reducing the 
overall normalization and therefore the implied mass de- 
position, M, also does not describe these observations. 
Clearly, the many aspects of these observations should 
offer some clue to the solution to this problem. 

3. Theoretical Challenges 

There are many ideas for possible mechanisms 
which would eithe r add or subtr a ct hea t from cooling 
flows (see e.g., iPeterson et all 1200 it iFabian et al.1 



iBohringer et all l2002t ISasaki fe Yamasakil 
20021 iMakishima et all 120011 IPeterson et all I: 
Mathews fe Brighenti [2003). .Promising scenari 



scenarios are 

document ed in these p rocee din gs and include dust 
coolin g l)Fabian et alJ 120011: iMathews fe Brighentil 
, cold cloud interface mixing (Begelman fe Fabiarl 
IFabian et all l200ll). interactions with cos- 
mic rays l|Gitti. Brunetti. fe Sett! 12002^ . heating 
through p rocesses associate d with relativistic AGN 
outflows jRosner fe Tuckerl Il989t iBinnev fe Tabor 
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esting explanations that might involve subtle errors in 
the application of Equation 1, such as a gross error in 
the plasma codes, are probably not relevant. This is 
likely not to be the case since some cataclysmic variables 
are observed to have X-ray spectra almost exactly lik e 
that in the top panel in Figure 2 ijMukai et al J 12003). 
This is understood by the idea that the accretion of 
material from the secondary to the white dwarf where 
it must cool down to the photospheric temperature 
essentially proceeds the same way as a classic cooling 
flow. 

Generally, there are three major requirements for any 
successful model. First, the average heating or additional 
cooling power should be comparable to the X-ray lumi- 
nosity from cooling flows. Most calculations have shown 
that this could be true for a number of different models. 
The second and somewhat more confusing requirement 
is that the mechanism should work more effectively at 
lower temperatures or should begin to operate right be- 
fore complete cooling occurs. The isobaric cooling time is 
proportional approximately to the second power of tem- 
perature. So if the cooling starts to be disrupted charac- 
teristically at ^ of the maximum temperature, why does 
it stop when it is | of the way to completion? Finally, 
a successful model must also add or subtract the heat at 
the right time at all spatial positions. The observed spa- 
tial distribution of the X-ray emission, however, is some- 
what complicated as we discuss in the following section. 
The global problem, however, is quite clear- the X-rays 
from the lowest temperatures are not present in the en- 
tire cooling flow volume. 

4. Comparison with other X-ray observations 

Some discussion on the nature of fitting spatially- 
resolved collisional X-ray spectra is warranted, since the 
comparison of these observations with those derived from 
other X-ray instruments can be complicated and was un- 
der frequent discussion at this conference. If the ICM 
is optically-thin to its own radiation and is in collisional 
equilibrium then its differential luminosity distribution 
is described by two functions, 
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where the first function fa describes the differential lu- 
minosity distribution as a function of temperature, T, at 
each solid angle element dil and the second function fa 
specifies an abundance map for each element. The mea- 
surement of other quantities, such as a density or vol- 
ume, are always indirectly derived from these two func- 
tions. Note also that the function, ^ t ^q , is not actually 
observable in practice. This is because a plasma con- 
sisting of multiple temperatures will emit line emission 
that can be separately measured from each temperature 
region but the continuum emission from each tempera- 
ture region will just contribute to one overall continuum. 
Even if the plasma were divided into different temper- 
atures each with their own metallicity, only one overall 
metallicity can be measured for any one spectrum. 

The ability to measure a temperature distribution, 4f , 
is determined by the number of counts in the spectrum, 
the resolution of the instrument, and fundamental atomic 
physics issues. For a high resolution spectrum with a 
large number of counts, the X-ray temperature distribu- 
tion probably cannot be divided into much more than 8 
or so temperature bins. This is because there are only 
so many strong emitting X-ray ions. The Fe K and L 
series, for example, have 10, but their overlap in specific 
temperatures ranges makes disentangling the tempera- 
ture distribution in great detail difficult. Line ratios be- 
tween emission lines from the same ion are only weakly 
temperature sensitive for the normal temperature ranges 
that the ion is produced. 

With a typical RGS cluster observation of about 30,000 
source photons, we have typically found that deviding the 
X-ray temperature range into 4 logarithmic bins usually 
provides a robust solution. With more photons, we ex- 
pect that a slightly more detailed analyses can be per- 
formed. With the non-dispersive CCD instruments of 
EPIC and Chandra less information can be extracted. 
Individual X-ray ions are not detected precisely, but their 
emission is blended significantly. It appears that usually 
determining much more than two temperatures (or two 
temperature ranges) is extremely difficult without a very 
large number of counts. 

The RGS observations place a different set of con- 
straints on the spatial distribution of the luminosity dis- 
tribution, ^ than EPIC or Chandra observations. The 
ability to do this, of course, depends on the spatial res- 
olution of the instrument provided enough X-ray pho- 
tons are collected. Therefore, Chandra can constrain the 
angular luminosity distribution at the arcsecond scale, 
whereas EPIC can constrain the distribution at the sev- 
eral arcsecond scale. RGS observations can utilize the 
cross-dispersion information (the distance photons are 
detected perpendicular to the dispersion axis of the spec- 
trometer) to measure the spatial distribution almost as 
good as the EPIC resolution. The shapes of line profiles 



at least in principle can constrain the distribution in the 
dispersion direction. 

The ability to measure the total two-dimensional func- 
tion, dTdn IS obviously complicated and probably is dif- 
ferent for every observation of every cluster. We may 
have to guess at what the true distribution or what 
the shape of this function is for cooling flows. Mak- 
ing conclusive statements about isothcrmality or a pre- 
ferred formulation of the temperature distribution is very 
challenging and depends heavily on the capabilities of 
that instrument. All of the instruments, however, have 
generally shown an average temperature decline. The 
RGS observations presented here show evidence for a 
steep differential luminosity distribution of order T 1 or 
T 2 if averaged over the entire cooling flow. The deep- 
est EPIC observations show steep multiphase distribu- 
tions (possibly power laws with indices greater than 
T 2 ) with overall temperature cont rasts that are prob- 
ably greater than a factor of two ijKaastra et alJ 120031 
iLewis et al.ll2003fl . Chandra thermal images have shown 
non-s pherical temperature va riations at levels less than 
50 % i|Sanders fe Fabianll2003j) . We might speculate that 
all the observations have given us a picture that the com- 
bination of an overall temperature decline and a steep 
multiphase distribution (T 3 ) at each radius combines to 
give an overall differential distribution near to what has 
been measured with the RGS (T 1 or T 2 ). This, however, 
has not been fully explored nor is it clear how much dif- 
ference there is between any of the clusters. 

The abundance map, is only beginning to be mea- 
sured since it depends heavily on getting the tempera- 
ture distribution correct. Furthermore, deviations from 
collisional equilbrium have yet to be fully explored. In 
addition, the optically-thin assumption as well could be 
violated by resonant scattering but th e effect is either 
observationally controversial (see e.g., iBohringer et ail 
2002] IXu et a,Ul2002l IChurazov etlu EoO?,) or could be 
different for every cluster due to varying velocity fields. 

5. Future Work 

As we have argued in the previous section, there is 
much more work to be done in clarifying the observa- 
tional situation in X-rays. Similarly, the connection of 
these observations to other wavelengths is very incom- 
plete. In particular, maps that correlate the quantity of 
missing X-rays with the detections of cold material wi ll 
be important l|Edge et all l2002t lEdee fe Fraverll2003|) . 
Theoretical studies should continue to build towards a 
final solution. 

It has been previously discussed that cooling flows of 
some form are necessary to cool galaxies to form stars 
before they get incor porated into much larger structures 
ijWhite fe Reeslll978|) . Unless the solution to the soft X- 
ray cooling flow problem is endemic to cluster scales, it 
could have far-reaching consequences for structure for- 
mation. 

JRP would like to thank the conference organizers and 
the attendees for a very successful conference. Work on 
the RGS at Columbia University and U. C. Berkeley is 
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